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Supplementary Material 

Consider an infinite, haploid population divided spatially into two denies, E x and E y . Each individual 
in the population is defined by two biallelic loci: a major locus A/ a which controls the phenotype and 
thus the fitness of the individual, and a modifier locus M/m which controls the switching rate between 
phenotypes A and a. Switching occurs only at the phenotypic locus, at a rate that is assumed to be the 
same in both directions. Selection acts only on the phenotypes A and a. The modifier locus M/m is 
assumed to be selectively neutral. 

Within each deme, the environment varies temporally between two states, Ti and T 2 . At each time- 
step, individuals first experience selection, followed by switching, recombination, and finally migration. 
The recursions representing the change in frequency at each generation are presented below. 

Equations describing the change in frequency at every generation. 

Denote by x\, x 2 , X3, X4 the frequencies of MA, Ma, mA, and ma in deme E x and let y\, y 2l V3, and 2/4 
be the analogous frequencies in deme E y . 

The fitness of phenotype i in temporal state j within deme I is denoted by Wiji (i in {A, a}, j in {1, 2}, 
I in {x, y}). The mean fitness in deme / is denoted by wi, I £ {x, y}. After selection, there is a process of 
switching between phenotypic states. There are two possible switching rates: ^m, associated with allele 
M and \i m associated with allele m. After selection and switching, the equations for the frequencies of 
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the four genotypes in the two denies are 



xi = (1 - mm) _ xi + hm^-x 2 
w x w x 

X 2 = MM - Xi + (1 - — X 2 

w x w x 

X3 = (1 - flraK^^ + Hm^-Xi 

w x w x 

~ A ix / \ o, ix / -t \ 

XA=^m _' X 3 + {l-fJ, m ) _' X 4 (1) 

W x W x 

n ,W A ,ky W aky 

2/1 = (1 ~Mm) _ 2/1 +Mm^!/2 

Wy Wy 

WA.ky , \ w a,ky 

2/2 = Mm— — 2/1 + (1 - MmH— » 

2/3 = (1 - Mm) _ " 2/3 + Mm^W 



WA,fej/ . /-. \ w a,ky 
2/4 = Mm" 2/3 + (1 - Mm)— 2/4, 



where the % and Wj, are the mean fitnesses in denies E x and -E y , respectively. (For example, w x = 
wajxXi + w a j x x 2 + wajxXz + w a .j X x 4 .) The indices j and k denote the current temporal state within 
each deme (j, k G {1, 2}). The index j refers to demc E Xl while the index k refers to deme E y . 

After recombination and migration, the frequencies x\ and y\ (i <E {1, 2, 3, 4}) at the start of the next 
generation can be written in terms of the post-switching frequencies Xi and j/i as 

x' x = {l-v) (Si - rOj) + v (j/! - rD 5 ) 
a/ 2 = (1 - i/) (x 2 + rflj) + v (y 2 + rDy) 
x' 3 = (1 - f) (£ 3 + rDj) + ^ (2/3 + ^-Dy) 

a/ 4 = (1 - v) (x 4 - rD x ) + v (y 4 - rDy) (2) 
y[ = v (x 1 - rD x ) + (l-u) (2/1 - rDy) 
2/ 2 = 1/ (i 2 + rD x ) + (1 - v) (2/2 + rL> 5 ) 
2/3 = f (&3 + rife) + 0--V) (2/3 + r^) 
y 4 = v (x 4 - rD x ) + (1 - v) (y 4 - rDy) , 



where D x and Dy are the linkage disequilibria in demc E x and E y , respectively (D x = x\x 4 — x 2 x 3 and 
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Dy = mm - mm)- 



Supplementary Figure Legends 
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Figure SI. Invasion trials showing modifier frequency as a function of generation time. The 

selection coefficients are all 0.1. Recombination and migration rates are 0.1. The environment switches 
deterministically in both demes every 10 generations. The stable switching rate in this case is 0.0271 
and is the selected switching rate for the modifier allele m. The switching frequency determined by the 
wild-type allele M is shown by dot color. We chose wild-type rates both above and below the modifier 
rate. In all examples, the modifier spreads to fixation in both demes. 

Figure S2. Robustness of the result to different symmetric selection pressures. All selection 
coefficients are equal and as presented in the legend. Recombination rate is 0.1. The expected time 
before an environmental change in both demes is 10 generations, and the environmental variability tp 
(the variance of the gamma distribution divided by the square of its mean) is 0.1. The plotted curves 
represent a fit to the data using a generalized additive model with penalized cubic regression splines 
Figure S3. Robustness of the symmetric selection results to different environmental mean 
waiting times. The selection coefficients are all 0.1. Recombination rate is 0.1. The environmental 
variability ip is 0.1. Panel A: The expected time before an environmental change is the same in both 
demes and presented in the legend. Panel B: The expected time before an environmental change is 10 
generations in deme E x . The expected time before an environmental change in deme E x is represented 
by line color. The plotted curves represent a fit to the data using a generalized additive model with 
penalized cubic regression splines. 

Figure S4. The stable switching rate as a function of migration when selection strengths 
differ between demes. The selection coefficients arc s\ = S2 = 0.4 and S3 = S4 = 0.1. Recombination 
rate is 0.1. The expected time before and environmental change in both demes is 10, and the environ- 
mental variability ip is represented by line color. The plotted curves represent a fit to the data using a 
generalized additive model with penalized cubic regression splines. 

Figure S5. Robustness of the asymmetric selection results to different environmental mean 
waiting times. The selection coefficients are Si — s 3 = 0.4 and s 2 = S4 = 0.1. Recombination rate 
is 0.1. The environmental variability ip is 0.1. Panel A: The expected time before an environmental 
change is the same in both demes and presented in the legend. Panel B: The expected time before an 
environmental change is 10 generations in deme E x . The expected time before an environmental change 
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in deme E x is represented by line color. The plotted curves represent a fit to the data using a generalized 
additive model with penalized cubic regression splines. 

Figure S6. Recombination and the stable switching rate. The selection coefficients are all 0.1 
and the environments switch with an expected period of 5 generations. The migration rate v = 0.05, 
and the environmental variability ip is represented by line color. The plotted curves represent a fit to the 
data using a generalized additive model with penalized cubic regression splines. 
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